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This paper summarizes the properties of strongly correlating liquids, i.e., liquids with strong 
correlations between virial and potential energy equilibrium fluctuations at constant volume. We 
proceed to focus on the experimental predictions for strongly correlating glass- forming liquids. These 
predictions include i) density scaling, ii) isochronal superposition, iii) that there is a single function 
from which all frequency-dependent viscoelastic response functions may be calculated, iv) that 
strongly correlating liquids are approximately single-parameter liquids with close to unity Prigogine- 
Defay ratio, and v) that the Active temperature initially decreases for an isobaric temperature up 
jump. The "isomorph filter", which allows one to test for universality of theories for the non- 
Arrhenius temperature dependence of the relaxation time, is also briefly discussed. 



I. INTRODUCTION 

After the initial reports in early 2008 of the existence 
of a class of strongly correlating liquids [l] [2] , these liq- 
uids were described in four comprehensive publications 
that appeared later in 2008 and in 2009 in the Journal of 
Chemical Physics [3]-[6]. This paper briefly summarizes 
the properties and characteristics of strongly correlating 
liquids as detailed in Refs. 3-6 and present a number of 
new computer simulations. We list a number of experi- 
mental predictions for strongly correlating liquids, focus- 
ing on glass-forming liquids since this volume constitutes 
the proceedings of the Rome conference held in Septem- 
ber 2009 (6IDMRCS). The main message is that the class 
of strongly correlating liquids, which includes the van 
der Waals and metallic liquids, are simpler than liquids 
in general. This explains, for instance, the long known 
observation that hydrogen-bonded liquids have several 
peculiar properties. 



II. STRONG VIRIAL / POTENTIAL ENERGY 
CORRELATIONS IN LIQUIDS 

Consider a system of N particles in volume V at tem- 
perature T. The virial W is defined by writing the pres- 
sure p is a sum of the ideal gas term NksT/V and a term 
reflecting the interactions as follows 



pV = NksT + W . 



(1) 



If [/(ri, Tat) is the potential energy function, the virial, 
which has dimension of energy, is given [7HTT] by 

W{r,,...,rN) = -l/3^r,-Vr.f/(ri,...,rA,). (2) 

i 

Equation ([T]) describes thermodynamic averages, but it 
also applies for the instantaneous values if the virial is 



defined by Eq. ^ and the temperature is defined from 
the kinetic energy in the usual fashion [Tfill . 

If AU is the instantaneous potential energy minus its 
average and AW the same for the virial at any given 
state point, the WU correlation coefficient R is defined by 
(where sharp brackets denote equilibrium NVT ensemble 
averages) 



R 



{AWAU) 



^{{AWr){{AUr)' 



(3) 



By the Cauchy-Schwarz inequality the correlation coeffi- 
cient obeys —1<R<1. We define strongly correlating 
liquids by the condition R > 0.9 [3]. The correlation coef- 
ficient is state-point dependent, but for all of the several 
liquids we studied by simulation [H HI [12] R is either 
above 0.9 in a large part of the state diagram, or not at 
all. 

Figure [l] shows two examples of constant-volume ther- 
mal equilibrium fluctuations of virial and potential en- 
ergy for two model systems, the standard Lennard- Jones 
(LJ) liquid and the Wahnstrom binary Lennard- Jones 
mixture [13]. In both cases there are strong virial / po- 
tential energy correlations. In (b) one sees striking dips 
in the potential energy; these dips reflect the existence 
of transient clusters in the liquid characterized by the 
same short range order as the crystal [T4|. During the 
dips virial and potential energy also correlate strongly. 
Actually, the correlation even survives crystallization [4] . 
Thus the property of strong virial / potential energy cor- 
relations is quite robust; even complex systems like bio- 
logical membranes may exhibit strong correlations [T5]. 

One way to illuminate the correlations is to plot in- 
stantaneous values of virial and potential energy versus 
one another in so-called scatter plots. Figure |2|a) shows 
an example of this with data taken from a simulation of 
the Kob- Andersen binary Lennard- Jones (KABLJ) liquid 
[16] . This has become the standard liquid for studying 
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FIG. 1: (a) Instantaneous normalized equilibrium fluctua- 
tions of virial and potential energy in the standard single- 
component Lennard- Jones liquid at constant volume (NVT 
simulation). W{t) and U{t) correlate strongly, (b) The same 
for the supercooled Wahnstrom binary Lennard- Jones mix- 
ture [13]; here W{t) and U{t) correlate strongly even during 
the formation of a so-called Frank-Kaspers cluster 1141. 



viscous liquid dynamics, because it is difficult to crystal- 
lize (this requires simulating for more than 100 microsec- 
onds (Argon units) [17 ). The "slope" 7 of the scatter plot 
gives the proportionality constant of the fluctuations ac- 
cording to 



AV^(t) ^ 7A[/(t) . 



(4) 



The number 7, which varies slightly with state point, is 
roughly 6 for the standard LJ liquid, roughly 5 for the 
KABLJ liquid, and roughly 8 for the OTP model studied 
below in Fig. [?[ 

Since viscous liquid dynamics consists of long-time vi- 
brations around potential energy minima - the so-called 
inherent states [18 - followed by rapid transitions be- 
tween the inherent states [I9l[20], it is interesting to study 




FIG. 2: (a) Virial / potential energy correlations for the Kob- 
Andersen binary Lennard- Jones liquid (1000 particles studied 
by Monte Carlo simulation, p — 1.264, T — 1.24 in standard 
L J units) . (b) Inherent state energies and virials of the simu- 
lation in (a); the correlation is still high. The slope 7 defined 
by AVK(t) = 7A[/(t) is slightly different, but comparable to 
that of the true dynamics. 



the inherent dynamics analogue of Fig. [2|a). This is 
done in Fig. [2|b), which gives the same simulation data 
after minimizing the configurations' potential energy us- 
ing the conjugate gradient method. The correlations are 
still present and the "slope" 7 doesn't change very much 
- even though the virial decreased by more than 60% go- 
ing from (a) to (b). This confirms the robustness of virial 
/ potential energy correlations. 

A convenient way to get an overview of a liquid's WU 
thermal equilibrium fluctuations at constant volume is 
to collect scatter plots for several state points in a com- 
mon diagram. Figure [s] (top) shows such a plot for the 
standard LJ liquid. Each state point is represented by 
one color. As in Fig. |2]the strong correlation is reflected 
in the fact that the ovals are highly elongated. For each 
value of the density the ovals form almost straight lines 
with slope close to 6 (in Ref. ^ it was shown that during 
and after constant-volume crystallization the system's 
scatter plots fall on the extension of the line). The bot- 
tom three figures show the correlation coefficient R (Eq. 



3 




Zero avg pres. 

p = 0.875 p = 0.840 P = 0.802 p = 0.761 

864 Lennard-Jones Particles 



p = 0.908^ 



-5500 -5000 -4500 -4000 

Instantaneous potential energy, U 



we reported simulations of 13 different model liquids. All 
liquids with van-der-Waals type interactions were found 
to be strongly correlating {R > 0.9), whereas models 
of the two hydrogen-bonding liquids water and methanol 
were not. Although much remains to be done by means of 
theory and simulation, it has now been established with- 
out reasonable doubt that liquids can be classified into 
two classes: (i) The class of strongly correlating liquids, 
which includes the van der Waals and metallic liquids; 
this liquid class has a number of regularities and simple 
properties, (b) All remaining liquids - the hydrogen- 
bonded, the covalently bonded, and (strongly) ionic liq- 
uids - which are much more complicated. 
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FIG. 3: (a) Scatter plot of the WIJ thermal equilibrium 
fluctuations at constant volume for the standard single- 
component LJ liquid, and (b) plots of various quantities as 
functions of temperature for the different densities studied. 
The full black line marks state points of zero average pres- 
sure. 



(|3])), the "slope" 7, and the average pressure as func- 
tions of temperature for the different densities. Clearly, 
both R and 7 are somewhat state-point dependent. At 
a given density R increases with temperature whereas 
7 decreases; at a given temperature R increase with in- 
creasing density. The thick black lines mark state points 
of zero average pressure. Note that the density effect of 
increasing R "wins" over the temperature effect of de- 
creasing R upon cooling at constant low pressure. Thus 
one expects higher correlations upon supercooling a liq- 
uid, which is an important observation when it comes to 
focusing on glass-forming liquids. 

How common are strong WU correlations? In Ref. 3 



III. CAUSE OF STRONG VIRIAL / 
POTENTIAL ENERGY CORRELATIONS 

Before discussing the consequences of strong virial / 
potential energy correlations we briefly reflect on the 
cause of the correlations. The starting point is the well- 
known fact [THTT] that for any liquid in which the parti- 
cles (of one or more types) interact with purely repulsive 
inverse power-law forces, v(t^ oc r~^, there is 100% cor- 
relation between W and \J\ W(i) = 7t/(t) where 



n 
3 



(5) 



From the values of 7 close to 6 observed for the LJ liquid 
one would expect that, if the LJ liquid somehow corre- 
sponds to an IPL liquid, the exponent n is close to 18. 
Although at flrst sight this may seem strange given the 
and T~^^ terms that enter into the deflnition of the 
LJ potential, a potential proportional to r~^^ does indeed 
give a good flt to the repulsive part of the LJ potential 
(Fig. [4] (a)). The reason that a much larger exponent 
than 12 is required is that the attractive term makes 
the LJ repulsion much steeper than that of the r~^^ term 
alone. Figure [4] (b) shows that both potential energy and 
virial fluctuations of the LJ liquid are well represented by 
those of an r~^^ IPL potential. 

In our flrst publications on strongly correlating liquids 
[TJ [2] it was suggested that the strong correlations de- 
rive from particle close encounters taking the intermolec- 
ular distance to values below the LJ potential minimum, 
at which the IPL potential is a good approximation. It 
quickly became clear, however, that this is not the full 
explanation; thus this can explain neither the existence 
of strong correlations in the crystal (above 99% at low 
temperatures [4 ), nor the existence of correlations at 
low pressures at which nearest-neighbor interparticle dis- 
tances fluctuate around the LJ potential's minimum dis- 
tance. Also, the original explanation is a single-pair ex- 
planation, which would imply that the strong correlations 
should be present as well in constant pressure ensembles. 
This contradicts our flnding that switching from constant 
volume to constant pressure reduces R from values above 
0.9 to values around 0.1 
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FIG. 4: (a) Approximation of the LJ potential by an effective 
inverse power law (IPL) potential oc r~^^. The blue dotted 
curve marks the IPL potential, which approximates the LJ 
potential well at small interparticle spacing. The red open 
circles mark the radial distribution function at a typical low- 
pressure state point. The difference between the LJ potential 
and the IPL potential is approximately linear in r; this fact 
forms the basis for the "extended inverse power law" (eIPL) 
approximation (Eq. (|6|) [4]|5j. (b) Two figures demonstrating 
that the LJ potential and its virial in their thermal equilib- 
rium fluctuations correlate strongly to the same quantities for 
the r~^^ IPL potential. 

References (4l and [5] detail the more complete expla- 
nation of the cause of strong correlations. The differ- 
ence between the IPL potential and the LJ potential 
is plotted in Fig. |4] (a) as the red dashed curve. The 
green dashed curve is a straight line, which approximates 
the red dashed curve well around the LJ minimum (i.e., 
over the entire first peak of the structure factor). Thus 
over the most important intermolecular distances an "ex- 
tended" inverse power law potential, eIPL, defined by 

^eiPL(r) = Ar-^^B^Cr (6) 

gives a good approximation to the LJ potential, vi,j{r) = 
'^eiPL(^)- It has been shown by simulation that the lin- 
ear "quark confining" term of the eIPL potential gives a 
contribution to the total potential energy that fluctuates 



little at constant volume [H [5] . Thus as regards fluctua- 
tions^ the pure IPL gives representative results. This ex- 
plains why the IPL approximation works so well and why 
the strong correlations disappear when going to constant 
pressure ensembles. This also explains why several IPL 
liquid properties are not shared by LJ-type liquids (e.g., 
the IPL equation of state is generally quite wrong and 
does not allow for low-pressure stable liquid states, and 
the IPL free energy and bulk modulus are quite wrong). 

While the eIPL approximation explanation of strong 
WU correlations for physically realistic cases, there are 
also strong correlations in the purely repulsive Weeks- 
Chandler- Andersen [21 version of the KABLJ liquid 
[ini [221 [23] . The slope 7 here varies quite a lot (from 5.0 
to 7.5) over the range of densities and temperatures in 
which 7 is fairly constant for the KABLJ liquid. Our sim- 
ulations show that the strong correlations for the WCA 
case is a single-particle-pair effect, not the cooperative 
effect that only applies at constant volume conditions, 
observed for LJ-type liquids. More work is needed to 
illuminate the correlation properties of this interesting 
(but physically unrealistic) potential. 



IV. ISOMORPHS: CURVES OF INVARIANCE 
IN THE PHASE DIAGRAM 

This section defines isomorphs and summarizes their 
invariants. As shown in Ref. 6 a liquid has isomorphs if 
and only if the liquid is strongly correlating. An isomorph 
is a curve in the phase diagram along which a number of 
properties are invariant. 

For any microscopic configuration (ri, ... ^r^) of a 
thermodynamic state point with density p, the "reduced" 
coordinates are defined by = p^^^r^. State points (1) 
and (2) with temperatures Ti and T2 and densities pi 
and p2 are said to be isomorphic [6 if, whenever two mi- 
croscopic configurations {r^i\ ... ,1*^'') and (rp\ ... ,r^^) 
have identical reduced coordinates, to a good approxima- 
tion they have proportional configurational NVT Boltz- 
mann probability factors: 

The constant C12 here depends only on the state points 
(1) and (2), not on the microscopic configurations. Iso- 
morphic curves in the state diagram are defined as curves 
for which any two state points are isomorphic. The prop- 
erty of having isomorphs is generally approximate - only 
IPL liquids have exact isomorphs. For this reason Eq. 
should be understood as obeyed to a good approximation 
for the physically relevant configurations, i.e., those that 
do not have negligible canonical probabilities [6 . 

Figure [s] illustrates Eq. ^ by checking the logarithm 
of this equation, where (a) gives simulation data for 
the KABLJ liquid. We consider a number of configu- 
rations of the state point with density and temperature 
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FIG. 5: Direct check of the isomorph condition for the KABLJ 
hquid (8000 particles), which is strongly correlating (a), and 
for the SPC water model (5120 molecules), which is not (b). 
For both liquids the consistency of the isomorph condition is 
checked by jumping from one to a different density and back. 
This works well for the KABLJ liquid but not for SPC water; 
details are given in the text. 



(pi,Ti) = (1.258,0.628) in standard LJ units. For these 
configurations the total potential energy was evaluated. 
In order to investigate whether the state point has an 
isomorphic state point at density p2 — 1.228 we scaled 
the simulated configurations of the first state point to 
density p2- For the scaled configurations the potential 
energies are plotted against the original energies of state 
point 1 (top figure). According to the isomorph defini- 
tion Eq. ([7|) the best fit slope gives the ratio between the 
temperatures of the isomorphic state points; in this way 
we estimate that T2 = 0.555. 

The right panel of Fig. (sja) investigates the consis- 
tency of this procedure by reversing it in order to check 
whether the original temperature T\ is arrived at. In- 
deed, when this is done one does find the original tem- 
perature to be 0.628. Two things should be noted. The 
first is the very strong correlation between scaled con- 
figurations, as required for having good isomorphs. The 
second notable fact is that the best fit lines do not pass 
through (0,0). This shows that the constant C\2 of Eq. 
([7|) is not unity, as it would be for an IPL liquid {C\2 is 
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FIG. 6: (a) A A particle radial distribution function of the 
KABLJ liquid for two isomorphic state points (left) and for 
the same temperature at the two densities (right). The iso- 
morphic state points have the same radial distribution func- 
tions, (b) The same for the A A incoherent intermediate scat- 
tering function at the k- vector corresponding to the first peak 
of the radial distribution function. The two isomorphic state 
points have the same dynamics (in reduced units, as used 
here) . 



determined by the contribution to the partition function 
coming from the linear term of the eIPL, and thus C12 
reflects the deviation from true IPL behavior). 

Figure |5|b) makes this "direct isomorph test" for the 
non-strongly correlating liquid SPC water, starting from 
temperature T\ = 200 K. From the slope of the left panel 
we find T2 = 179.6 K. When the reverse jump is per- 
formed, however, one does not come back to the initial 
state point, but to a predicted temperature of 166.36 K. 
This shows that water does not have isomorphs, consis- 
tent with the fact that it is not a strongly correlating 
liquid. 
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For the practical identification of an isomorph in the 
phase diagram the above method may be used. Alterna- 
tively, it has been shown [6 that to a good approximation 
isomorphs are characterized by 



Const. 



(8) 



Here 7 is the above discussed "slope" characterized by 
AW{t) = jAU{t). As shown in Ref. 6 this quantity may 
be calculated to a good approximation from equilibrium 
fluctuations via the expression (giving the least-squared 
linear-regression best-flt slope of WU scatter plots, com- 
pare Appendix B of Ref. |3]) 
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(AWAU) 

am') ' 



(9) 



Several physical quantities are invariant along a 
strongly correlating liquid's isomorphs to a good approx- 
imation. These include: 1) Thermodynamic properties 
like the excess entropy (i.e., in excess of the ideal gas en- 
tropy at same density and temperature) and the excess 
isochoric speciflc heat, 2) static averages like radial dis- 
tribution function(s) in reduced coordinates, 3) dynamic 
quantities like the reduced diffusion constant, viscosity, 
and heat conductivity, time-autocorrelation functions in 
properly reduced units, average reduced relaxation times, 
etc. 

Figure [6] shows results of simulations of the KABLJ liq- 
uid at two isomorphic state points (left subflgures) and 
isothermal state points (right subflgures) of the AA par- 
ticle radial distribution functions and the AA incoherent 
intermediate scattering functions, respectively. These flg- 
ures conflrm the prediction that isomorphic state points 
have identical static distribution functions and identical 
dynamics. 

Since the isochoric speciflc heat is an isomorph invari- 
ant, this quantity should be a function of p^/T for a 
strongly correlating liquid. Figure [7| conflrms this for 
the Lewis- Wahnstrom OTP model consisting of three LJ 
spheres [24 . 

The theory further predicts that jumps between two 
isomorphic state points should take the system instanta- 
neously to equilibrium, because the Boltzmann statisti- 
cal factors of two isomorphic state points by deflnition 
are proportional [6.. More generally, isomorphic state 
points are equivalent during any aging scheme. We re- 
cently showed that the isomorph concept can be used to 
throw light on the concept of an effective temperature 
|26. . In particular, the theory implies that for strongly 
correlating liquids the effective temperature after a jump 
to a new (low) temperature and a new density, depends 
only on the new density (Fig. (sja)). We showed also 
that this does not apply for the non-strongly correlat- 
ing monatomic Lennard- Jones Gaussian liquid, conflrm- 
ing the general conjecture that strongly correlating liq- 
uids have simpler physics than liquids in general. Fig- 
ure [s] (b) shows a result conflrming the flnding of Ref. 
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FIG. 7: cv per particle for the Lewis- Wahnstrom model [23] of 
ortho-terphenyl consisting of three LJ spheres arranged with 
fixed bond lengths. The left panel shows the raw simulation 
data, The right panel shows the same data replotted as func- 
tion of p^'^ /T in which the exponent 7.9 was determined from 
the proportionality between equilibrium virial and potential 
energy fluctuations [12]. This model follows the Rosenfeld- 
Tarazona prediction of cy varying with temperature as T~^/^ 



26 that the effective temperature concept for a strongly 
correlating liquid makes good sense physically. On the 
X-axis the inherent state energies of given state points 
are shown as the system fell out of equilibrium. The ar- 
rested phase is characterized by an effective temperature 
Teff, which can be calculated in the standard way from 
the violation of the fluctuation-dissipation theorem [26] . 
On the y-axis is shown the inherent energies found from 
an equilibrium simulation with temperature equal to Teff 
for the corresponding arrested phase. The red points 
give data for the strongly correlating KABLJ liquid, the 
green points give data for the non-strongly correlating 
monatomic Lennard- Jones Gaussian liquid. Clearly, the 
latter system fell out of equilibrium by freezing into a 
part of phase space that is not characterized by Teff. 



V. THE EQUATION OF STATE OF A 
STRONGLY CORRELATING LIQUID 

This section shows that the Helmholtz free energy for 
any strongly correlating liquid is of the form 



Fex/iV = T^{Tf{p))+g{p). 



(10) 
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FIG. 8: (a) Virial versus potential energy for the KABLJ 
liquid during a temperature quench and a crunch 26 . A 
crunch increases the density and keeps the temperature con- 
stant; this is equivalent to first an instantaneous jump along 
an isomorph (green curve) to the right density followed by a 
temperature quench (black curve). The inset shows that the 
crunch and the quench have the same fluctuation-dissipation 
violation factors, i.e., result in the same effective tempera- 
ture, (b) Inherent state energies for several state points of 
the KABLJ liquid as the system falls out of equilibrium and 
freezes versus the equilibrium inherent state energy at the 
effective temperature (details are given in Ref. ^26 ) . 



For simplicity we shall not indicate "excess" quantities 
explicitly, V is used as variable instead of the p, and the 
is ignored since it is fixed; thus we shall prove that 



Since Cy = T{dS/dT)v = T(j)'-^{x){dx / dT)v , if we de- 
fine m{x) = 02(^)/^i(^)5 T{dx/dT)y = m{x). 
Along an isochore this implies that dx/m{x) = dlnT. 
When integrated along the isochore this gives h{x) = 
j^^dx' /m{x') = InT + a(y), implying ex.p{h{x)) = 
Tf{V) where f{V) = exp(a(V)). In other words, lines 
of constant x have constant Tf{V). Since x is merely 
used for labeling the isomorphs, this means that we 
can redefine X as follows X = Tf(V). Integrating now 
S = (j)i{x) = —{dF/dT)y along an isochore gives 

F = - [ Mx')dT'+g{V) = -- r Mx')dx'+g{V). 

JTo ^ Jxo 

(13) 

The integral is some function of x, and we have thus 
derived Eq. ( pT| . 

As mentioned the slope 7 may vary with state point. 
The equation of state gives information about which vari- 
ables 7 may depend on: 

S = - (f^) = -HTfiV)) - T^P'{Tf{V))f{V) 

(14) 

implies that 

U = F + TS = -T^i;'{Tf{V))f{V)+g{V). (15) 
We also have 

y = - {w)^ = -T^^'iTf{V))f{V)-9'{V). 

(16) 

Combining these equations leads to 

w = (u - aiv)) - va'iv) . (i7) 

Summarizing, 

W = -f{V)U^h{V), (18) 

where 



F = Ti^{Tf{V))+g{V). 



(11) 



Suppose a given strongly correlating liquid's isomorphs 
are labeled by the variable x, i.e., that its isomorphs are 
curves of constant x for which x = x{T,V). Isomorphs 
are curves of constant (excess) entropy, as well as curves 
of constant (excess) Cy [6 . This means that for some 
functions (/)i(x) and (/)2(^) one may write 



S = (/)i(x) , Cv • 



(12) 



7(F) = 



d\nf 
dlnV 



and 



HV) = -giV) 



dHfg) 
dlnV 



(19) 



(20) 



Thus if 7 varies with state point, it may only depend on 
volume [6]. As shown in Ref. 6 this result is consistent 
with the original experimentally based formulation of the 
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so-called "density scaling" due to Alba-Simionesco and 
co-workers [27] . Figure [sj^b) shows that 7 is not rigor- 
ously constant along isochores as predicted by Eq. (18), 



although it varies only of order 10% when temperature 
is tripled. This serves to emphasize that isomorphs are 
only approximate constructs and so are their predicted 
invariants (only IPL liuqids have exact isomorphs). 
The equation of state Eq. (18) is of the Mie-Griineisen 



form for the excess variables (excess pressure: W/V^ and 
excess energy: U) ^ 



(21) 



In the Mie-Griineisen equation of state for solids ^{V) = 
—dlnuj/dln V where a; is a phonon mode eigenfrequency, 
U is the vibrational potential energy, and tp{V) relates 
to the volume derivative of the energy per atom (i.e., 
of the energy of the force-free configuration about which 
the vibrational motion occurs). Further discussion of the 
relation of strong WU correlations to the Mie-Griineisen 
equation of state is given in Ref. ^Si 



VI. SOME EXPERIMENTAL PREDICTIONS 
FOR STRONGLY CORRELATING 
GLASS-FORMING LIQUIDS 

Strongly correlating liquids have a number of since long 
described properties. For instance, since the melting line 
in the phase diagram is an isomorph 0, strongly corre- 
lating liquids have several invariants along their melting 
lines, including the radial distribution function and di- 
mensionless transport coefficients. Such regularities have 
been observed in simulation and experiment; we refer to 
the reader to Ref. 6 for more details. This section focuses 
on predictions for highly viscous liquids, for which the 
strong-correlation property implies several experimental 
predictions. 

1. Density scaling 

In the last decade, in particular since 2005, many 
papers appeared dealing with the so-called density 
scaling, which is the finding that for several glass- 
forming liquids the relaxation time r at varying 
pressure and temperature is some function of the 
quantity p^/T, in which the exponent 7 is an em- 
pirical fitting parameter: 



F(pyT) . 



(22) 



Neither the function F nor 7 are universal. The 
isomorph theory [6l [12] shows that all strongly cor- 
relating glass formers obey density scaling with the 
exponent 7 given by the equilibrium fiuctuations 
at one state point (Eq. ([9|) (provided 7 is fairly 
constant over the relevant part of phase space). 



This has not yet been tested experimentally, but it 
is consistent with the finding that density scaling 
works well for van der Waals liquids [29 , but not for 
hydrogen-bonded liquids [ 30H33) . Meanwhile, den- 
sity scaling has been shown to apply in computer 
simulations of strongly correlating liquids with 7 
given by Eq. ([9| to a good approximation [T2j |34] . 

2. Isochronal superposition 

Isochronal scaling is the further, fairly recent find- 
ing that for varying pressures and temperatures 
the dielectric loss as a function of frequency de- 
pends only on the loss-peak frequency [35l [36] , 
This is trivial if the liquid obeys time-temperature- 
pressure superposition (TTPS) in which case noth- 
ing changes. But many liquid do not obey TTPS, 
and for such liquids isochronal superposition is a 
new and striking regularity that works generally 
for van der Waals liquids, but rarely for hydrogen- 
bonding liquids [36] . Since both the relaxation time 
and the entire relaxation spectrum are isomorph 
invariants, isochronal superposition must apply for 
any strongly correlating liquid: If temperature and 
pressure for two state points are such that their re- 
laxation times are the same, the two points must 
belong to the same isomorph and thus have same 
relaxation time spectra for any observables, in par- 
ticular the dielectric loss as function of frequency 
should be the same. 

3. Frequency- dependent viscoelastic response func- 
tions 

There are eight fundamental complex, frequency- 
dependent linear thermoviscoelastic response func- 
tions like, e.g., the frequency-dependent isochoric 
or isobaric specific heat, the frequency-dependent 
isobaric expansion coefficient, and the frequency- 
dependent adiabatic or isothermal compressibility 
[37 . Standard linear irreversible thermodynamic 
arguments, where the Onsager relations play the 
role of the Maxwell relations of usual thermody- 
namics, show that there are only three independent 
frequency-dependent response functions. If more- 
over stochastic dynamics is assumed as is realis- 
tic for highly viscous liquids [38 , there are only 
two independent response functions [39MT] . For 
strongly correlating liquids the further simplifica- 
tion appears that there is just a single indepen- 
dent response function [TJ |37l [42] . Since there are 
explicit expressions linking the different response 
functions (depending on the ensemble considered 
[37 ), this can be tested experimentally. Unfortu- 
nately it is difficult to measure thermoviscoelastic 
functions properly; to the best of our knowledge 
there are yet no reliable data for a complete set of 
three or more such response functions on any liquid. 

4. The Prigogine-Defay ratio: Strongly correlating liq- 
uids as approximate single-parameter liquids 
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After many years of little interest the Prigogine- 
Defay (PD) ratio [43ti45] has recently again come 
into focus in the scientific discussion about glass- 
forming liquids [37l [46ti49] . From a theoretical 
point of view the PD ratio is poorly defined since 
it involves extrapolations from the liquid and glass 
phases to a common temperature [ST, "42] . It is pos- 
sible to overcome this problem by modifying the PD 
ratio by referring exclusively to linear response ex- 
periments; here the traditional difference between 
liquid and glass responses is replaced by a difference 
between low- and high-frequency values of the rel- 
evant frequency-dependent thermoviscoelastic re- 
sponse function [37j. In this formulation, the prop- 
erty of strong virial / potential energy correlations 
manifests itself as a PD ratio close to unity. Actu- 
ally, an extensive compilation of data showed that 
van der Waals bonded liquids and polymers have 
PD ratios close to unity [50 - even though as men- 
tioned the traditional PD ratio is poorly defined, 
there are good reasons to assume that it approxi- 
mates the rigorously defined "linear" PD ratio. 

The theoretical developments of Refs. [21 011371 show 
that in any reasonable sense of the old concept 
"single-order-parameter liquid", strongly correlat- 
ing liquids are precisely the single-order-parameter 
liquids. The isomorph concept makes this even 
more clear: State points along an isomorph have so 
many properties in common that they are identi- 
cal from many viewpoints. In the two-dimensional 
phase diagram this leaves just one parameter to 
classify which isomorph the state point is on; thus 
a liquid with (good) isomorphs is an (approximate) 
single-parameter liquid. Note that this is consis- 
tent with the old viewpoint that single-parameter 
liquids should have unity PD ratio [43-45 . 

5. Cause of the relaxation timers non-Arrhenius tem- 
perature dependence: The isomorph filter 

Since the relaxation time r is an isomorph invariant 
for any strongly correlating liquid, any universally 
valid theory predicting r to depend on some phys- 
ical quantity must give r as a function of another 
isomorph invariant (we do not distinguish between 
the relaxation time and the reduced relaxation time 
since their temperature dependencies are virtually 
identical). This gives rise to an "isomorph filter" 
[6], showing that several well-known models can- 
not be generally valid. For instance, the entropy 
model cannot apply in the form usually used by 
experimentalists: r ex exp(C/5'conf^) where C is a 
constant and 5'conf is the configurational entropy; 
it can only be correct if C varies with density as 
C (X . Likewise, the free volume model does not 
survive the isomorph filter. If the characteristic 
volume Vc of the shoving model (predicting that 
r (X ex.-p{VcGoo/kBT) [5TH53]) varies with density 
as Vc cx this model is consistent with the iso- 



morph filter. 

6. Fictive temperature variations following a temper- 
ature jump 

Any jump from equilibrium at some density and 
temperature to another density and temperature 
proceeds as if the system first jumped along an iso- 
morph to equilibrium at the final density and then, 
immediately starting thereafter, jumped to the fi- 
nal temperature (Fig. (sja)): The first isomorphic 
jump takes the system instantaneously to equilib- 
rium. This property, which applies for all strongly 
correlating liquids, means that glass- forming van 
der Waals and metallic liquids are predicted to 
have simpler aging behavior than, e.g., covalently 
bonded liquids like ordinary oxide glasses [26 . 

In traditional glass science the concept of "fictive 
temperature" is used as a structural characteristic 
that by definition gives the temperature at which 
the structure would be in equilibrium [54H59] . For 
any aging experiment, in glass science one assumes 
that the fictive temperature adjusts itself mono- 
tonically from the initial temperature to the final 
temperature. Consider, however, a sudden tem- 
perature increase applied at ambient pressure. In 
this case there is first a rapid thermal expansion 
before any relaxation takes place. This "instan- 
taneous isomorph" takes the system initially to a 
state with canonical (Boltzmann) probability fac- 
tors corresponding to a /oiL'er temperature. In other 
words, immediately after the temperature up jump 
the system has a structure which is characteristic 
of a temperature that is lower than the initial tem- 
perature. With any reasonable definition of the fic- 
tive temperature, this quantity thus initially must 
decrease during an isobaric positive temperature 
jump - at least for all strongly correlating liquids. 

VII. SUMMARY 

The class of strongly correlating liquids includes 
the van der Waals and metallic liquids, but excludes 
the hydrogen-bonded, the covalently bonded, and the 
(strongly) ionic liquids. Due to their "hidden scale in- 
variance" - the fact that they inhering a number of 
IPL properties - strongly correlating liquids are sim- 
pler than liquids in general. Strongly correlating liquids 
have isomorphs, curves along which a number of physical 
properties are invariant when given in properly reduced 
units. In particular, for glass-forming liquids the prop- 
erty of strong virial / potential energy correlations in the 
equilibrium fiuctuations implies a number of experimen- 
tal predictions. Some of these, like density scaling and 
isochronal superposition, are well-established experimen- 
tal facts for van der Waals liquid and known not apply 
for hydrogen-bonded liquids. This is consistent with our 
predictions. Some of the predicted properties have not 
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yet been tested, for instance that the density scahng ex- 
ponent can be determined by measuring the hnear ther- 
moviscoelastic response functions at a single state point, 
or that jumps between isomorphic state points take the 
system instantaneously to equilibrium, no matter how 
long is the relaxation time of the liquid at the relevant 
state points. - We hope this paper may inspire to new 
experiments testing the new predictions. 
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